Non Poissonian statistics in a low density fluid 
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Our interest goes to the collisional statistics in an arbitrary interacting fluid. We show that even in the low density 
limit and contrary to naive expectation, the number of collisions experienced by a tagged particle in a given time 
does not obey Poisson law, and that conversely, the free flight time distribution is not a simple exponential. As 
an illustration, the hard sphere fluid case is worked out in detail. For this model, we quantify analytically those 
deviations and successfully compare our predictions against molecular dynamics simulations. 



The purpose of the present paper is to address, for an ar- 
bitrary low density fluid in a stationary state, basic properties 
that have been under-appreciated or overlooked, which bear 
upon the collisional statistics: What is the probability dis- 
tribution P{N , t) of the number of collisions JV suffered in 
equilibrium by a tagged particle over a given duration t? Con- 
versely, what is the probability distribution function of the 
free flight time, P^{t), again for a tagged particle? Such 
fundamental questions, relevant in their own right, have also 
consequences on the evaluation of transport coefficients and 
when transposed to the related realm of granular gases |[lll2l, 
directly quantify dissipation. Our message is that intuitive ex- 
pectation fails -except in highly untypical cases- as far as the 
collisional statistics is concerned, for an interacting fluid in 
and also, arguably less surprisingly, out of equilibrium. All 
results reported are new, together with the kinetic theory tech- 
niques used. 

Whereas in a dense fluid, velocity correlations and hydro- 
dynamic effects are responsible for a non trivial collisional 
statistics, one could naively expect that in the dilute limit 
where collisions become uncorrected and molecular chaos is 
enforced iH 0], collisional events define a Poisson process 
P{uj) so that P{Af,t) = exp{-Ljt){u;t)-^ /Ml, where uj is 
the mean collision frequency (i.e. {J\f) /t oj at long times, 
where the brackets denote an ensemble average). The corre- 
sponding free flight time distribution would then be 
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However, as we shall see below -and this seems to have been 
ignored in the literature lUt]- such a point of view is flawed. In 
essence, the collision frequency for a particle with velocity v 
depends on v (it generically increases with v = |v|), which in 
turn induces correlations between successive collision times. 
After general considerations that encompass equilibrium and 
non equilibrium steady states, we will show that the colli- 
sional statistics is generically non Poissonian. This is what 
prompted us to focus on the simplest -analytically tractable- 
interactions and consider the equilibrium hard sphere fluid for 
illustrative purposes. Such a model is one of the most useful 
paradigms in statistical mechanics and has played an essential 
role in the development of the theory of liquids |6'] . It is one 
of the simplest system exhibiting a phase transition. Remark- 
ably, it does not only provide a valuable theoretical starting 



point, but also enjoys direct experimental realizations It 
is therefore surprising that such a well studied system yields 
non trivial properties in a limit where little would have been 
expected. We will see that at late times, when Lot ^ M, the ef- 
fect of such correlations is that P{M, t) is of Poissonian form 
but with a renormalized frequency ti;/V2 instead of u. This 
result holds irrespective of space dimension. In addition, ex- 
plicit and accurate results will be reported for the number of 
collision cumulants {M^)c- Our analytical predictions will be 
compared to numerical simulations. 

For the sake of simplicity, we begin the analysis by the free 
flight time distribution P^^^ (t). The evolution of a tagged par- 
ticle in a large stationary homogeneous fluid defines, in the 
low density limit, a Markov process where the transition rates 
can be computed from a linearized Boltzmann equation, see 
e.g In simplified situations, the velocity dependent colli- 
sion rate r{v) can be computed analytically; the inset of Fig. [T] 
shows an iUustrative example for hard spheres in equilibrium. 
Velocities are expressed here in rescaled units, and for equilib- 
rium situations, the velocity distribution function reads, with 
d the space dimension 
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Out of equilibrium, is a stationary measure. In any case, the 
mean collision frequency oj follows from <j>: 



(v) r{v) dv = (r(w)) . 



(3) 



From the Markovian property, it follows that the conditional 
probability of having a free flight time t given a velocity v 
reads Pppj, (r|v) = r{v) exp(— r(z;)r). To proceed further and 
obtain P^ (t) from some average of P^ (^l v), attention must 
be paid to the fact that the relevant weight to use is not the 
velocity distribution 0(v) itself, but the velocity on collision, 
r{v)<f>{v)/uj. Theprefactorr(w) biases the distribution toward 
more energetic events and accounts for the fact that in a given 
time interval, a particle with a larger than typical velocity will 
experience more collisions. We therefore obtain 
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This expression explicitly differs from the result reported in 
Jgt], where r^/w is replaced by r (in other words, the weight 
used in is and not r<j)/Lu). To see why such an ap- 
proach is incorrect, one can compute the mean collision time 
(r) = / tP^^{t) dr, that should be equal to l/uj. This is in- 
deed the case with the distribution given in (HI, whereas the 
formula of Ref. |p| gives (r) = (l/r), which differs from 
1 / (r) = l/uj llOll . More importantly, upon neglecting the 
V dependence of the rate r (i.e. assuming r = oj), the inte- 
gral in (HJi is readily integrated and yields expression ([TJ for 
Ppp^. It has been shown that a w-independent collision rate 
corresponds to particles interacting via an inverse power law 
pair potential with exponent 2d — 2 |4(], which defines the so- 
called Maxwell model ifTlll . a particularly convenient frame- 
work in kinetic theory. Maxwell molecules are nevertheless 
highly untypical and for any other fluid, r depends on v so 
that (IDi cannot be exponential. We therefore conclude here 
that the collisional statistics is in general non Poissonian, ex- 
cept for Maxwell molecules where successive collisions turn 
out to be uncorrected. We will clarify below the conditions 
for the occurrence of correlations, and show that while the v 
dependence of r is a necessary condition for non Poissonian 
behavior, it is in general not sufficient. 

After the previous qualitative remarks, our goal is to quan- 
tify the deviations for Poissonian behavior, and to this end, we 
hereafter consider the prototypical hard sphere model where 
the frequency r{v) takes the form i^: 
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where iFi denotes a confluent hypergeometric function of the 
first kind. Although a closed-form expression cannot be ob- 
tained for Pppj (r) due to the lack of simplicity of the collision 
rate r{v), finding the large r behavior calls for a saddle point 
approximation for the integral appearing in Eq. dU, which 
yields 
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Interestingly, to leading order, r is distributed exponentially, 
as naively expected [see Eq. ([T]i], but with a renormalized 
rate L>jj\f2. The validity of expression (|6]l is illustrated in Fig. 
[T] which displays results of numerical simulations. We have 
checked that the molecular dynamics data in Fig. [Uprecisely 
coincide with the numerical integration of Eq. (|4]i for all ve- 
locities (not shown). 

We also note here that very similar considerations hold for 
the distribution of path lengths: for a particle with velocity v, 
the free flight distance (FED) t travelled in a time t\& I ~ vt, 
so that PppJ^lv) = v-^P,^iev-^\w). 
the counterpart of Eq. (|4]i 



Pppp then follows from 
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Figure 1 Free flight time distribution as a function of ujt, on a linear- 
log scale, for a two dimensional hard disc gas (d = 2). The circles 
correspond to the results of molecular dynamics simulations at den- 
sity p = O.OAa~^ where a is the discs diameter, with A*" = 1000 
particles. The dotted line shows the prediction of Eq. ^ while Eq. 
Q is shown with a dashed line. The inset shows the tagged colli- 
sion frequency r{v)/u) versus v, also for d = 2. The value 1 /\/2 at 
V = 0, indicated with an arrow, is d-independent. 



The resulting probability density is not a simple exponential, 
at variance with a claim sometimes found in the literature (see 
e g- lfl2i [Tsll ). However, a saddle point computation akin to 
that leading to ^ provides here the long distance behavior 
of Pppp, which is governed by the minimum of the function 
r{v)/v iH. This leads a) to P,,„i£) cx cxp(-a;(u)-i^/V2) 
for £ much larger than the mean free path, and b) to the re- 
mark that for the Maxwell model with a v independent rate 
r, the minimum of r{v)/v is reached for w — > oo and van- 
ishes, which leads to a non-exponential large £ Umiting be- 
havior in . 

We now turn to the related but more complex question of 
the number of collisions. We introduce the joint probability 
f{v,Af, t) of having velocity v and having suffered Af col- 
lisions in a time window t, for our tagged particle. The cor- 
responding time evolution, again in the dilute limit, follows 
from the linear Boltzmann-like equation 



9t/(vi,A/',t) = y dv2 y dCT(vi2 ■B)0{-vi2 ■ B) 

[/(vr,AA-l,t)<^(vr)-/(vi,AA,i)0(v2)] , (8) 

where 6 is the Heaviside function, vi2 is the relative veloc- 
ity, is a unit vector and the superscript refers to pre- 
collisional velocities: v" = vi — (vi2 • a-)a and Vj* = 
V2 + (vi2 • a-)a-. The latter equation encodes a full descrip- 
tion of the collisional statistics for the tagged particle, in the 
low density limit. Here again we stress that such an analytical 
approach can be extended to other interaction potentials, lead- 
ing to the already mentioned remark that, apart from Maxwell 
molecules, it does not admit a Poisson solution. Moreover, 
the above equation does not admit a stationary solution in the 
long time limit, due to the time dependent behavior of the col- 
lision number M. More precisely, we expect the large time 
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dependence to be exponential, as a consequence of the lin- 
ear character of the equation. For analytical progress, it turns 
convenient to introduce the generating function / through 



J\f=Q 
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Of course, upon summing (|8]l over all possible values of 
Af, in the equilibrium state one recovers X]^=o fi'^i-^) ^ 
f{v,0,t) = 0(v), the Maxwell-Boltzmann distribution 
It then appears that the cumulant generating function /i(A), 
such that in the large time limit 



t(-ir ^ 
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is the largest eigenvalue of an evolution operator that straight- 
forwardly follows from ^ ifTsll . Furthermore, /i(A) is directly 
related to the large time behavior of P{J\f, t) through its large 
deviation function tt, defined as 

P{M,t) ""^^ e*'^^") (11) 
where n = J\f/t. Indeed, tt is the Legendre transform of yu: 

7r(n) = min(/i(A) + Xn) . (12) 

The quantity /i(A) therefore bears an important physical in- 
formation, and has been the technical focus of our study. 

From perturbation theory, we have obtained the behavior of 
/i at large A in the form 



l) + 0(e-2^) , 



which implies that 
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For this Poissonian behavior 7'((jj/\/2), the large deviation 
function easily follows: 



7r(ri) = n — nlog(n\/2/w) — lo/\/2 . 



(15) 



We note that Eq. ( fT4] l with A/^ = is compatible with the time 
integral of the leading exponential order of free flight time 
distribution given in as it should. 

However, the dependence embodied in Eqs. (fTsT l and (fl4] i 
follows from a large A expansion and only holds for M tot 
(hence n/oj ^ 0). It is therefore not indicative of the typi- 
cal behavior, for which it is more useful to consider the low 
A limit. We then dwell on the remark made after ^ that for 
A = 0, we have /(v, 0, t) — 0(v), which leads to the ap- 
proximation scheme where /(v, A, t) is taken as a Gaussian 
with temperature T{X) and the requirement that T(0) = 1 
iflill . This approximation is the lowest order of a more sys- 
tematic expansion (see ifTsIl ). but it provides a useful and rea- 
sonably accurate information given its simplicity. The idea is 
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Figure 2 Large deviation function n{n) characteristic of the long 
time behavior of the probability P(N ,t) of suffering N collisions 
in a time t (and defined by ^{M /t) ~ t"^ log P{J\f,t) for cot » 
1). The symbols correspond to molecular dynamics measures for 
a hard disc system with A'^ = 1000 particles at two different but 
low densities. The solid line shows the Gaussian result following 
from Eqs. l |16l > while the dashed line is the large deviation function 
7r(n) = n — nlog(n/aj) — uj associated to the Poisson law Viuo). 
On the graph, the time scale has been set by the choice uo = 1. 



to project the Boltzmann equation ([8]) governing the evolution 
of /(v, A, t) onto the first two velocity moments. This leads 
to 
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It can be checked that /i(oo) = —uj/\/2, as implied by Eq. 
(fT3] l. The corresponding first three cumulants follow from 
([Tol l and read, irrespective of dimension d 
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These values are compared to molecular dynamics simulation 
data in TableJ] Moreover, the large deviation function 7r(n) of 
collisions, which follows from ( fT2b . appears to be in excellent 
agreement with its molecular dynamics counterpart, see Fig. 

El 

At this point, it proves instructive to consider the Lorentz 
gas with only one mobile particle and a collection of spherical 
fixed scatterers. One readily gets r{v) cx v, but the collisional 
statistics is nevertheless Poissonian. The reason is that the 
velocity modulus of the mobile particle is constant along the 
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TABLE I Cumulants for the number of collisions M from molecular 
dynamics simulations (performed on a two dimensional system with 
reduced density pa^ — 0.04), and comparison with both the Poisson 
'P{<jj) result and Gaussian approximation 
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trajectory: in other words, there is no thermalization. It there- 
fore appears that non Poissonian behavior arises from two key 
properties, that induce colHsional correlations: first, the colli- 
sion rate depends on the velocities, and, second, the particle 
thermalizes to some stationary non singular measure. These 
criteria also apply out of equilibrium. 

In conclusion, we have shown that the statistics of the num- 
ber of collisions J\f experienced by a tagged particle in a low 
density homogeneous and stationary fluid, in or out of equi- 
librium, is a subtler quantity than it might seem. Our general 
statements have been illustrated with the hard sphere fluid. 
For the distributions of both J\f and the related free flight 
time -for which several incorrect results may be found in the 
literature-, we have quantified the corresponding non Pois- 
sonian behavior which follows from the simple physical in- 
gredient that a particle with a high velocity statistically col- 
lides more often than a typical particle. A key quantity in 
the theoretical analysis is the cumulant generating function 
/i(A), which can be computed explicitly for large A and ap- 
proximately for small A. The Gaussian ansatz worked out 
here can be considered as the lowest order of a systematic ex- 
pansion. The resulting analytical predictions have been con- 
fronted against molecular dynamics numerical simulations, 
with a very good agreement. These numerical results show 
that the deviations from Poisson behavior V{uj), although 
not dramatic -which may be the reason why they are under- 
documented in the literature- are nevertheless clearly observ- 
able. In particular, we have obtained the a priori surprising 
result that for long times, the distribution of J\f is Poisso- 
nian, but with a "dressed" rate w/\/2. Conversely, the dis- 
tribution of free flight time r is exponential, with a behavior 
(X exp(— ijjt/V2). 
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